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ABSTRACT 

Preliminary experiments in the numerical prediction of large scale 
sea-surface temperature anomalies are made using a 10-level primitive 
equation model with 300-km horizontal resolution covering a rectangular 
basin in the North Pacific. The model is first integrated over an 11- 
year period to statistical equilibrium using time dependent wind and 
thermal forcing. The monthly normal climatology, generated by the 
model in this long term integration, is then used along with observed 
sea-surface temperature anomalies to define the initial state. 

The results show that inclusion of salinity in the model may pro- 
duce more accurate predictions of high latitude cold anomalies. The 
most accurate of the predictions removed all of the anomalous diffusion 
terms in the governing equations. Data for the observed temperature 
anomalies were not available at depth, thus initialization of the 
anomalies below the surface was superficial. Consequently, the effects 
of vertical advection were not fully realized, indicating a need for 
future experimentation with a more accurately defined initial profile 
of the anomalies in the initial conditions. 
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I. INTRODUCTION 



For today's meteorologist to make accurate, long range forecasts the 
understanding of the dynamics of atmospheric behavior is imperative. In 
the recent meteorological past, the importance of air-sea interactions 
has been realized and has given the forecaster another potential tool with 
which to improve his predictions. A better understanding of the inter- 
actions between the oceans and the atmosphere will also assist in ocean 
acoustic forecasting and other Naval operations. 

Like the atmosphere, the oceans experience seasonal, monthly, and 
daily variations in temperature and circulation patterns. Further investi- 
gation of these variations reveals the existence of anomalous behavior 
within the mean temperature and circulation patterns. It is believed 
that these variations influence directly the atmospheric weather patterns 
through the exchange of moisture, momentum, and heat across the air-sea 
interface. This thesis is concerned with the prediction of these 
anomalous oceanic temperatures and currents. 

During the past two decades it has become apparent through the studies 
of J. Namias and the late J. Bjerknes that the atmospheric response to 
sea-surface temperature anomalies may be worldwide. Bjerknes (1966, 1969, 
1972) has shown that large-scale, sea-surface temperature anomalies exist 
in the equatorial region and that these anomalies produce a response in 
the atmosphere which is not limited to that region alone. Likewise, 

Namias (1959, 1969, 1972b) discusses the possibility that such a response 
is also present in the middle latitudes. Furthermore, Namias (1969) 
observed that during several years when intense sea-surface temperature 
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anomalies were present in the North Pacific Ocean , there were distinct 
weather patterns over the United States. Because of the above pioneer- 
ing findings of Namias and Bjerknes, the North Pacific Experiment 
(NORPAX) was established in which a cooperative scientific effort is 
being made to understand the causes of these anomalies and their rela- 
tionship to the atmosphere. 

Numerical integrations of coupled ocean-atmosphere models have been 
reported by Manabe (1969) , Bryan (1969) , and Wetherald and Manabe (1972) . 
By comparing the results of their coupled model with those of the ocean 
and atmosphere models integrated separately, these authors showed the 
significant role ocean circulation plays in determining the state of 
the atmosphere. 

In turning our attention to the oceans, the first attempts at large 
scale ocean predictions was based upon the use of a prediction equation 
for the sea-surface temperature (Namias, 1965). In Namias* method, 
temperature T is a function of x,y,z, and time t. The local rate of 
change of temperature can be represented by 



3T 3t 3t 3t 

3T = - u 37- V 37- W 3l +q 



(1) 



where q is the rate of change of temperature from non-advective pro- 
cesses. Now if we let ( ) denote the long term normal state and (*) the 
anomaly, then the monthly mean can be expressed by 

T=t+T* ,u=u+u* , v = v + v* ,w=w+w* 

and 

q = q + q' . 
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Then Eq. (1) takes the form: 



d_ 

3t 



(T + T’) 



- (u+u ’) 



3 (T+T') 
cTx 



(v+v ' ) 



3 (T+T ' ) 

3y 



- (w+w' ) 



3 (T+T 1 ) 
3z 



+ q + q’ 



( 2 ) 



Assuming the long term normal temperature satisfies the equation 



3T 



- 3t - 3 t - 3t 



^ = - u -k v-r w — + q 

3t 3x 3y 3z 



(3) 



subtracting (3) from (2) gives the equation for the temperature anomaly: 



3T ' 

3t 



- ,, 3t' - 3t' - 3t; 

[(u+u') — +(v+v') — + (w+w') — ] 
3x 3y 3 z 



. , 3T, , 3T a , 3t. , 

[u' -r— + v' ^— + w' -T— ] + q' . 
3x 3y 3z 



(4) 



Namias ' first predictions (1965) retained only the advection of the 

3 T 3 T 

normal temperature by the anomalous currents, (u 1 + v* 7 ^), in Eq, 

(4). The anomalous currents (u*, v') were obtained from the anomalous 



atmospheric wind stress observed during the period of the "forecast" 
and the long-term normal temperatures (T) from observations collected 
by the U. S. Weather Bureau. Since this prediction used data actually 
observed during the forecast period, it is not a prediction in the true 
sense and is referred to as a "semiprediction". Although this simple 
method produced results which correlated positively with the observed 
temperature changes, Arthur (1966) demonstrated that additional terms 
should be included in the prediction equation. As a result, Jacob (1967) 
performed a semiprediction neglecting only the vertical advection terms 
and the anomalous heating. A quantitative comparison of his results 
indicated an improvement in the prediction. 
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More recently, Namias (1972a) has shown that the sea-surface tempera- 
ture patterns in the North Pacific can be predicted with some success 
using normal surface currents to advect the anomalies (neglecting all but 
-u ^ , — v on the right side of Eq. (4)). The results of his study 

(Figs. 1, 2, 3 , and 4), which will be discussed in detail later, show 
the importance of horizontal advection. These results plus the lack of 
temperature data at depth clearly indicates the need for a completely 
dynamical model. Therefore, it is the objective of this thesis to predict 
the evolution of the anomalous temperatures using a numerical model which 
contains all the advection terms in (4) and a parameterization of mixed 
layer dynamics as well. 

To make a comparative analysis, the period covered by ‘this prediction 
will be the same as that used by Namias (1972a) . Sea-surface temperature 
anomaly data were also obtained from Namias 1 20-year normal observations. 
It should be noted, however, that this normal data is from the years 1947 
to 1966 and is not necessarily representative of the prediction period. 
Furthermore, due to the aforementioned lack of data at depth, we assume 
an anomalous temperature profile which decreases linearly from its value 
at the surface to zero at 150 meters. NORPAX scientists have shown that 
anomalies exist at depths greater than 150 meters; however, this depth 
was chosen to correspond more closely to Namias* study which assumed the 
temperature anomaly to be confined strictly to the surface. Further 
studies could be conducted varying the vertical profile of the tempera- 
ture anomaly to determine the importance of the initial condition. The 
initial anomalous current is geostrophic and its vertical average is 
zero. Unlike Namias* method which employed the use of observed 
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climatological data during the entire forecast period , this method uses 
model-generated climatology only to define the initial conditions. Con- 
sequently, following initialization the anomalous and mean fields are 
completely free to dynamically interact. 
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II. THE NUMERICAL MODEL 



A. BACKGROUND 



The basic model used in the following predictions is that of Haney 
(1974) as modified by Davies (1975) . "The model is based upon the hydro- 
static and Boussinesq approximations. This means that the ocean is 
assumed incompressible, and the density is replaced by a constant every- 
where except in the hydrostatic equation where it is multiplied by the 
acceleration of gravity. A major simplification is the neglect of 
salinity and the assumption that the density is a linear function of 
temperature alone."^ 

The equations which govern Haney’s model are given by 



-1 



du _ 

dt p^acoscj) 9A 



■I? +fv 



u . _ 9 2 u 

7 * v sin <f) + A V u + K , 

acos<J> m j 2 



(5) 



dv _ -1 9p 

dt p a 9<J> U 
o 

9p 

■ ‘ pg ' 



acos<}> 



2 8 2 v 

u sin <p + A V v + k 



m 



9z 



2 ' 



9w 1 r 9u 9 , j. , , 

■ 5 — + T [-TT- + ~T- (v COS <}>) ] =0 , 

9z acostp 9A 9<J> 

dT „ „2 9 2 T , /m . 

— = A V T + K - + 6 (T) , 

dt H ^ 2 c 

9z 



P = p [1 - a (T - T ) ] , 

o o 



where the representation of the symbols is as follows: 

t time 
X longitude 
(J> latitude 



(6) 

(7) 

(8) 

(9) 

( 10 ) 



Haney, R. L. , "A Numerical Study of the Response of an Idealized 
Ocean to Large Scale Heat and Momentum Flux," J. of Physical Oceanog ., 
V. 4, No. 2, p. 147, 1974. 
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z height 

ft angular speed of the earth's rotation 
f coriolis parameter 
a radius of the earth 
g acceleration of gravity 
u eastward velocity component 
v northward velocity component 
w vertical velocity component 
T temperature 

T constant reference temperature 
p density 

p density of water at reference temperature 
a coefficient of thermal expansion 



B. BOUNDARY CONDITIONS 

After Haney (1971) and Davies (1975) , the determination of the upper 
boundary condition in the thermal energy equation is computed by calcu- 
lating separately the incoming solar radiation (Q^) , the net upward long- 
wave radiation (Q ) , the latent heat exchange (Q ) , and the sensible heat 
B E 

exchange (Q ) . Together, these fluxes define the net downward heat flux 
H 

as represented by 



pressure at 10 meters, C the fraction of cloud cover, and 0 the 
Stefan-Boltzmann constant. The exchange of latent heat and sensible heat 
across the ocean-atmosphere interface were computed from 




(ID 



The back radiation Q is calculated from 

B 



Q„ (T ) = . 985 ( . 39 - .05e 1/2 ) (1 - .6C 2 ) C T 4 
Bo a o 



( 12 ) 



where T q is the sea surface temperature, e^ the atmospheric vapor 




(13) 
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and 



Q = p C V C (T - T ) , 

*H K a D 1 1 p o a 



respectively , where p was the density of the air, C the drag co- 

efficient, |v| the wind speed at anemometer level, L the latent heat 

of vaporization, the atmospheric pressure, e^ (T q ) the saturation 

vapor pressure at the sea-surface temperature, the specific heat 

of air at constant pressure, and T the atmospheric temperature at 

a 

10 meters. Incoming solar radiation (Q^) was calculated with 



Q = S (.74 - .60 
I o 



(15) 



where S was the solar insolation incident on a horizontal surface at 
o 

the top of the atmosphere. 

The quantities S Q ' ^ ' e a ' T a anc ^ l V i are P rescr ik ec * functions 

of latitude (v depends on longitude — see below) and time taken from 

atmospheric climatological data. The fluxes Q , Q , and Q are 

BE H 

additionally dependent on the computed sea-surface temperature, T 

o 

Therefore, an anomalous sea-surface temperature produces a corresponding 
anomalous heat flux which has the tendency to remove the temperature 
anomaly by virtue of the opposing sign in Eq. (11) . 

"In addition to the thermal forcing, the ocean circulation was also 
forced by surface winds given on a time-dependent basis. Assuming an 
average geostrophic inflow angle of 10° and a frictional reduction of 
velocity by 10%, one could derive the surface stress from the surface 
geostrophic winds. If u and v are the eastward and northward surface 
winds, respectively, after inflow and frictional influences are considered, 
then 



V = 



, 2 

(u + v ) 



16 



and the wind stress components are 




Therefore, the wind stresses are completely specified since u(y,t) 
and v(x,t) are completely independent of the ocean temperatures and 
currents in this experiment. Since the quantities which define the 
atmospheric forcing are derived from climatological data, this thesis 
does not concern itself with the response of the ocean to anomalous 
atmospheric forcing. 

Lastly, the ocean bottom is assumed to be flat which requires the 
vertical velocity to vanish there. Thus integrating the continuity 
equation up to the surface and applying a free surface boundary condi- 
tion, one could obtain a prediction for the height of the surface. 

However, the present method assumes that the ocean surface is not a free 
surface but rather a balanced surface at which w = 0. Therefore, since 
w = 0 at the top and bottom, the vertical average of (8) shows that the 
vertical average current is non-divergent and consequently can be de- 
scribed by a streamfunction. The streamf unction is predicted from the 
vorticity equation which is derived from the vertical average of the 
equations of motion. 

C. RESOLUTION AND DATA 

To improve the resolution, the 6-level model employed by Haney (1974) 
was expanded to include 10 levels (Davies, 1975). The vertically expanded 



2 

Davies, R. W. , A Numerical Parameterization of Wind Mixing in a 
Time-Dependent Baroclinic Oceanic General Circulation Model , M. S. Thesis, 
Naval Postgraduate School, Monterey, California, 1975, pp 56. 
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model now has fixed levels at 10, 30, 60, 100, 150, 225, 350, 700, 1500 



and 3000 meters. 

As mentioned previously, to compare the results with those of Namias 
(1972a) , sea-surface temperature anomaly data covering the same period 
were obtained from J. Namias of Scripps Institution of Oceanography. 

The surface temperature data thus obtained covered only the area 20°N to 
60°N, 130°E to 110°W in the North Pacific and the grid spacing intervals 
for the latitude and longitude were 5° and 5° respectively. Therefore, 
Namias 1 data were expanded to the equator and to 65°N by extrapolating 
the anomalous data to zero at the most distant boundaries. Additionally, 
the data over the entire domain were interpolated to fit the 2.0° lati- 
tude and the 2.8° longitude grid spacing of the present model. 

As a result of the expansion, the predictions outside the original 
data field should be ignored. Consequently, for the convenience of the 
reader, an overlay (Fig. 5) with the locations of ocean vessel stations 
November (N) , Papa (P) , and Victor (V) as well as the site of the NORPAX 
"Pole” experiment, has been superimposed on the normal fields and each 
related prediction. The superficial land mass boundaries of the NORPAX 
grid are represented by the shaded area and the original data field by 
the area between the dashed lines. 

D. INITIAL CONDITIONS 

To completely define the initial state of the ocean for dynamical pre- 
dictions, a considerable amount of model climatology was used along with 
the available observed data. The model climatology was generated by in- 
tegrating the model over an 11-year period starting from a stratified 
initial state and using atmospheric data which contained the observed 
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normal seasonal data. This model climatology was used for initializa- 
tion and verification as described below. To provide continuity in 
the atmosphere wind and thermal forcing, the annual mean and amplitude 
and phase of the first two harmonics were calculated from atmospheric 
climatological data as defined by Davies (1975) . The atmospheric data 
were averaged along each grid latitudinal line, except for v, the north- 
ward component of the surface geostrophic wind, which was averaged 
along each meridian. 

The resulting initial normal currents and the 30, 60, 90 day normal 
currents at level 2 (30 m) are shown in Figs. 6, 8, 10, and 12, where 
"initial 11 is defined as October 15th of year 11. As seen in Figs. 6 and 
12, there is an almost negligible change in the normal currents during 
that period. The corresponding initial, 30, 60, and 90 day normal tem- 
perature fields are shown in Figs. 7, 9, 11, and 13. Although there is 
little change in the temperature gradient in the middle latitudes, the 
southward movement of the main baroclinic zone indicates a cooling trend 
which was expected during the winter months. 

To begin the predictions, the total field must be defined in the 
initial conditions. Thus, the normal part of the initial conditions is 
obtained from the model-generated climatology while the anomalous part 
is obtained from observations. The observed sea-surface temperature 
anomalies are extrapolated linearly from observed values at the surface 
to zero at level 5 (150 meters) . This rather arbitrary procedure of 
defining the anomalous temperatures below the surface is done because 
of the lack of sub-surface temperature data. 
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The total current (normal plus anomaly) can be represented by 



W = W + W 

where (^) is the vertical average current 



( 17 ) 



i r * 

, [~ J ( )dz] , ( ) is the 

-H 



vertical shear current, [( )-( )], and H is 
the total ocean depth. Expanding (17) into its normal and anomalous 
parts one obtains 



W = W + W* = W + \V ,+ W + W' 



(18) 



where W and \V are the normal part of the vertical average current 
and the normal part of the vertical shear current respectively. The 

A 

anomalous part of the vertical average current (W 1 ) is assumed to be 

* 

zero and the anomalous part of the vertical shear current (\V 1 ) is the 
geostrophic current calculated from the above temperature anomaly using 
the thermal wind relation. Since 



W = — — Ik x Vp , 
of 
o 



g P~ f 



taking the partial derivative with respect to z produces 



a\v 



9z p f 
o 



Ik x V (-pq) 



Substituting (10) for p 



p =-p a T + constant , 
o 



where the gradient of p becomes 



Vp = -p a g V T . 
o 
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Therefore, 



fk x VT (thermal wind relation) . 

9 z f 

It is now assumed that the anomalous part of the vertical shear current 
satisfies this relationship, 



9z 



f 



Ik x VT' • 



Then integrating from -H to some depth z 

/ z 

T'dC (19) 

-H 

* * 
where \V* (-h) is unknown. If the vertical average of \V* is taken, 

then this quantity can be set equal to zero by definition. Thus, from 

(19) 



o z 



* 

w* 



= | f \v'(- H)dz IkxV i / (/ T'dOdz = 0 , 



-H 



-H -H 



or 



o z 



\V 1 ( — H ) + 



^ lk x V | f ( ^ T'dO 



dz = 0 



( 20 ) 



* 

Then (20) can be solved for W'(-H) and substituting into (19) to give 
the result 



\Y 1 (z) 



og 

f 



Ik x V 



z 

( J T'dz 

-H 



o z 




( 



/ 



-H 



T'd£) dz] 



( 21 ) 



In Eq • (21) f is not allowed to become smaller than its value at 

13.25°N. This geostrophic approximation to the vertical shear current 
anomaly is then used as the initial condition in the present model. 

Having defined the initial temperature and currents by combining 
the model generated normal fields with the observed anomaly fields , the 
model is run for 90 days of simulated time. During this time, the total 
temperatures and currents interact and evolve together. At any given 
time later (e.g., 60 days) the predicted anomaly is defined to be the 
total field predicted for that time minus the model generated normal 
for that same time. 



E. THE PREDICTION EQUATION 

Since several cases were studied using variations of the thermodyna- 
mic prediction Eq. (9) which governs the total temperature field, a 
detailed analysis of that equation follows. 

Neglecting the convective adjustment (6^(T)) in the analysis, (9) 
can be written 



dT 

dt 




T + 




( 22 ) 



Expanding (22) 



9T 9t 9t 9T 

9t + U 9x +V 9y +W 9z 



2 2 2 

9 t 9 i 1 ^ 9 T 

A( — + — ) + K — 

8x 8y 8z 



(23) 



22 



and substituting further for the normal and anomalous parts , (23) takes 



the form 



3t 3t* - 3t - 3t* , 3t . 3t' - 3t - 3t* 

3l + 3¥ = ' [u 3^ + u 3^ + u 3^ + u 3^ +v ^T + v 97 



, 3t , 3t’ - 3t - 3t' , 3t , 3t'. 

+ v' - 77 - + v' 7 — + w 7 — + w 7 — + w' 7 — + w' -r— ] 

dy dy dz dz dz dz 



+ M i!| + 3fr' + + lc[ iff + sfl.', . 



3x 2 3x 2 3y 2 3y 2 



2 2 
3z 3z 



(24) 



Setting T', u ! , v* , and w 1 equal to zero, then the normal field can be 
represented by 



3t r - 3 t - 3t . - 3t, . , r 3 2 T . 3 2 t, . , 3 2 t 
3z^ = - u 3x + V 3y" + W 1Tz + M 7T + TT 1 + K 7~2 * 



(25) 



3 x 3 y 



3 z 



Subtracting (25) from (24), the equation governing the anomalous tempera- 
tures becomes, 



3t' 

3t 



- 3 t' ,3t ,3t' - 3t* ,3t , 3 t' 

- [u - 7 T- + U -7T- + U Sr- + V - 77 — + v 'rr- + v 'yr— + W 



3 x 



3 x 



3 x 



3y 



3y 



3y 



3T* 

3z 



— 3 1 

. 3t t 3T'] + a[-2— 

+ wSr- + w'yr- \ 2 

3 z 3 z 3 x 



2 2 
3 TS 3 T* 

+ ] + K r 



3y 



3 z 



(26) 



The following investigation will reveal that initially \y ,# VT' = 0 
by analysis of (21) . If we let "D" represent the depth at which the tem- 
perature anomaly becomes zero as defined earlier, then the temperature 
anomaly (T") at some level z < 0 is 



T 1 (x ,y , 0) - T (x , y , -D) 
D 



T' (x , y , 0) - T' 
- z 
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where T* (x , y , 0) is the sea-surface temperature anomaly and T’(x,y,-D) 



is the temperature anomaly at D = 150 meters. Since T* (x,y,-D) 
definition, then 



T' = T' (x,y, 0) (1 + J) . 



Substituting (27) into term (a) of (21) gives 



/ 



T 



dz 



-H 



F (z) T 1 (x,y,0) , 



where 



F(z) 



z 

/ <1 * i)d C 

-H 



and (^) is a dummy variable representing the vertical coordinate, 
wise, (28) can be substituted into term (b) of (21) to produce 



o z 

i / ( / T ' d « az 

-H -H 



i I° fm 

-H 



T* (x,y,0)dz 



A 

= F T‘ (x,y, 0) . 



Therefore, (21) becomes 

* Otq ^ 

W' = -f Ik x V [F(z) T' (x,y , 0) - FT' (x,y,0)] 

or by factoring T' (x,y,0) 

* CtQ a 

\V = J Ik x V h (T' (x,y, 0) (F (z) - F) ] . 
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0 by 



(27) 



(28) 



Like- 



(29) 



From this, it is clear that F(z) and F are not functions of x and y and 
can be brought outside the "del" operator allowing 



\V« = ^ (F (z) 



F) Ik x V h T' (x,y,0) 



(30) 



Consequently, 



\V'*V T' = ^-(F(z)-F) [ Ik x V T' (x,y , 0) *V T' (x,y,0) (1 + ^r) ] 

Hr H H u 



= ^-(F(z)-F) (1+ |) [ Ik x V h T' (x,y,0)»V H T' (x,y,0)) 



(31) 



and since term (a) of (31) is parallel to the anomalous temperature 
field and term (b) of (31) is perpendicular to the same field, then 

* /*v 

W 1 • V T ' = 0 o Since \v' =0 , we have 
H 

W' • V h T' = 0 , 

9t * 3t 1 

which allows u an< ^ v< ’3^’ to ne< ?l ectec ^ i n (26) . Figures 14 
and 15 clearly show that initially the anomalous currents are geostrophic 
and parallel to the anomalous isotherms. It should be noted that these 
currents need not remain parallel to the anomalous isotherms or geos- 
trophic during the prediction period; however , careful examination of 
the predictions indicates they do. It is not surprising that this type 
of motion remains geostrophic, but it is surprising that it remains equi- 
valent barotrophic. Therefore, for the purpose of this interpretation, 

we will assume that W 1 • V T 1 is negligible during the forecast period. 

H 
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3t 1 

Additionally , w'tj— can be further evaluated. From the continuity 



equation (8) , 



' 3w 
3z 



+ 



w 



0 , 



then 



3w 

3z 



+ V H 




+ V H 



w 



0 



where 



3w 

3z 



+ 



W = 0 , 



implying 



3w 1 
3z 



+ V 



W' 



0 . 



Thus integrating up from the level z to the surface z = 0 , 



w 1 




w* dz + w 1 (0) 



(32) 



where w 1 (0) = 0 as defined by the boundary conditions. Therefore, (32) 
can be written 



= f V 

J H 



* W' dz , 



(33) 



and since V • \V 1 = 0 , V • \y' = V • \y ' allowing the substitution of 



(30) into (33) as 



w' = ag 



/ V H 



[--' V r - F - Ik x V u T' (x,y,0) ] dz 
r H 



(34) 
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Since the coriolis parameter (f) is a function of y, and T' (x,y,0) is 



not a function of z, Eq. (34) becomes 



w' = V H # X V H T'(x,y,0) F (z) ] 



(35) 



where 

o 

F (z) = f (F(z) - F) dz . 
z 

Hence, evaluating (35) 



v - 3T '^.' y,0) i ^ j) ? (z)) 



9y 



9x 



(36) 



If we denote the geos trophic current, integrated from the level z to the 
surface z = 0 by (Ug, Vg) , 



Ug(z) 



qg 9 t* ( x,y , 0) 
f 9y 



and 



Vg (z) = 



qg 9T" (x,y , 0) - 
_ — F 



then (36) can be written 



w 1 



“»> 




(Vg)] . 



Therefore, 



w 1 



ccg 9f 9 t’ (x,y, 0) ~ 
^2 9y 9x 



6 

f 



Vg (z) 



and so. 



w' 



9t ' 

9z 



6 

f 



Vg (z) 



9t ' 

9z * 



(37) 



(38) 
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The approximation for the vertical velocity w 1 applies after integration 
begins as well as initially. This result shows that w' is the vertical 
motion produced by the divergence of the goestrophic current. Thus, a 
northward moving layer of surface water (Vg > 0) converges to produce 
sinking motion below the surface (w 1 < 0) , while southward moving surface 
water (Vg < 0) diverges and produces rising motion below the surface 
(w 1 > 0). As a result of the above evaluated terms of (26), the final 
prediction equation for the anomalous temperatures becomes 



3t' 

St 



r - 3 t' - 3t* , 3t j 

- [u + V *r— + U 1 *r — + 



3x 



3y 



3x 



4 5 



3T - 3t* 
it - + w 
3y 3z 



6 



6 

f 



Vg (z) 



3t 

3z 



- f Vg(z) ||'] + 



A[ 



a 2 T ' 

9x 2 



♦ ^§'l 

9y 2 



+ K 



3 2 t ' 

3z 2 



(39) 



7 



8 9 



The first two terms on the right-hand side of Eq. (39) represent 
the horizontal advection of the anomalous temperature by the normal 
current. These are the only terms used by Namias (1972a) to make his 
prediction. Terms 3 and 4 are the horizontal advection of the normal 
temperature by the anomalous current. Term 5 is the vertical advection 
of the anomalous temperature by the normal current. The vertical advec- 
tion of the normal temperature by the anomalous current and the vertical 
advection of the anomalous temperature by the anomalous current are 
terms 6 and 7 respectively. Term 8 is the horizontal eddy diffusion of 
the anomalous temperature and term 9 is the vertical eddy diffusion of 
the anomalous temperature which also includes the anomalous part of the 
surface heat flux (see discussion below (15)). 
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The motivation behind the inclusion of all the terms of (39) is, of 
course, to make a more accurate prediction,. As a result of Jacob's (1967) 
work, which showed an improvement over Namias' predictions (1972a) by 
including terms 3 and 4 of (39) , and since vertical advection is so im- 
portant in the theory of the permanent thermocline, the implication is 
that all the advection terms should be included when making a prediction 
of the sea-surface temperature anomaly. 
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III. RESULTS 



Before discussing the results of the predictions generated by the 
complete dynamical model, a brief description of the observed anomalous 
temperature fields and the results of Namias ' (1972a) predictions fol- 
lows. The initial, 30, 60, and 90 day observed anomalous temperature 
fields, as seen in Figs. 15, 16, 17 and 18, are averages over three 
months centered about the months of October, November, December, and 
January respectively. That is, the initial anomalous temperature field 
is the average of the anomalous temperatures for the months September, 
October, and November (i.e., the Fall season) of 1971. 

From Figs. 15 through 18 two cold anomalies and one warm anomaly 
appear to dominate the area of interest. The cold anomaly near the 
western boundary is observed to decrease in amplitude almost linearly 
during the first 60 days with a small southward movement. The double 
cell cold anomaly near the northeast boundary decreases in amplitude 
in the western part but maintains its relative amplitude and position in 
the eastern part. Finally, the warm anomaly near the center of Fig. 15, 
increases in amplitude with the leading edge of the maximum anomalous 
temperature contour moving eastward some 13° longitude. 

If one observes Namias' 90 day prediction of the warm anomaly (Fig. 
4), as described above, the easterly movement of the anomaly, as 
advected by the normal winter sea-surface currents (Fig. 3), is a good 
approximation. However, the amplitude of the anomaly in his prediction 
remains relatively constant, when it was observed to increase. In an 
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attempt to explain this discrepancy, Namias calculated the surface heat 
fluxes, but the results obtained indicated a negligible (and at times 
opposed) role of these fluxes. It is therefore suggested that the over- 
simplified use of the prediction equation, as described earlier, and the 
lack of data at depth is responsible for not predicting accurately the 
amplitude of the anomaly. 

Attempting to improve the predictions made by Namias (1972a) , the 
completely dynamical model of this thesis was run for a corresponding 
90-day period using the initial anomalous temperature field (Fig. 15) 
and the model generated climatological data along with the boundary con- 
ditions as described earlier. Four cases using variations in the 
boundary conditions and this prediction equation are discussed below. 

A brief description of the four cases will be given first, followed 
by a qualitative discussion of the results obtained. 

CASE 1: 

In CASE 1 the anomalous temperature prediction contained all the 
terms implied by Eq. (39). The results of CASE 1 are shown in Figs. 19- 
24. 



CASE 2: 

In CASE 2 the model calculated the 90 day prediction as in CASE 1 
except in the surface thermal boundary condition, Eqs. (12), (13), and 

(14) , the normal climatological sea-surface temperatures were substituted 
for the predicted sea-surface temperatures. This effectively removed the 
anomalous part of the surface heat flux. The results of CASE 2 are shown 
in Figs. 25-30. 
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CASE 3: 



CASE 3 is similar to CASE 2 except normal climatological temperatures 
were substituted in place of the predicted temperatures in the diffusion 

terms of (9) . This effectively eliminates the anomalous diffusion terms, 

2 2 2 
d T 1 9 T* 9 T 1 

[A (— — — + - — — ), K - ] , from (39), The results of CASE 3 are shown in 

9x 9y 9z 

Figs. 31-36. 



CASE 4: 

CASE 4 is similar to CASE 3 except normal climatological currents 



(u,v) were substituted for the predicted currents in the diffusion terms 
2 2 2 9 ^ 

(A V u , K — and A V v , K *) of the primitive Eqs . (5) and (6) . 

9z 9z 2 

This has the effect of removing the anomalous diffusion terms (A V u 1 / 

9 ^u * 2 9 ^v * 

K r and A V v', K — ) from the equations governing the anomalous 

9z 9z 

currents. The results of CASE 4 are shown in Figs. 37-42. 



Concentrating on the results of the anomalous temperature predictions 
of CASE 1 (Figs. 15, 20, 22, and 24), one observes that all the anomalies 
appear to decrease in amplitude with the cold anomalies disappearing com- 
pletely by the end of that 90 day prediction. The amplitude of the warm 
anomaly initially was found to be .6°C located at 168°W, 40°N and is 
reduced to .1°C after 90 days with its center at 162°W, 37°N. As a 
result of the above observations, cases 2, 3, and 4 were designed to 
determine the cause of the amplitude decrease or disappearance of these 
anomalies. Furthermore, the western edge of the warm anomaly is observed 
to move eastward while the eastern edge remains relatively stationary. 
Below 20°N some anomalous activity apj:>ears to have been generated; 
however, since this lies outside the area of initial observed data, it 
is disregarded and further discussion of this region will be omitted. 
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As in CASE 1, the cold anomaly at the northeast boundary disappears 
in CASE 2 (Figs. 15, 26, 28, and 30) . Analyzing the initial vertical 
profile of the total temperature field, it was found that the lapse rate 
was unstable thereby implying that convective adjustment (term 6^ (T) in 
(9) ) , is the cause of the cold anomaly vanishing. This rapid disappear- 
ance of the cold anomaly immediately raises the question of the impor- 
tance of salinity. If salinity were included in the present model, it 
is speculated that the cold anomaly could maintain itself. That is, the 
cold anomaly at the surface would not necessarily be unstable if a (nega- 
tive) salinity anomaly were present to stabilize the water in the column. 
Another possibility is that salinity is negligible, but the cold anomaly 
extends all the way to the ocean floor. However, this explanation seems 
less likely since it indicates the presence of such a large anomaly of 
heat content. 

Investigation of the initial vertical temperature profile of the 
total field in the vicinity of the cold anomaly near the western boundary 
reveals that it too is convectively unstable. However, this anomaly damps 
less rapidly in CASE 2 than does the cold anomaly at the northeast 
boundary which would indicate its decrease in amplitude is not entirely 
due to convective adjustment. Additionally, the amplitude of the central 
warm anomaly was found to decrease less (to a value of .2°C) than in 
CASE 1. Except for the above two differences, the 60 and 90 day predic- 
tions of cases 1 and 2 are relatively similar. Therefore, with respect 
to the observations of the first two cases, it is concluded that the 
anomalous surface heat flux induced by the anomaly itself is not the 
major cause of the damping of the anomalies, and that convective adjust- 
ment is the major cause of the cold anomalies vanishing. The effects of 
diffusion will be discussed in cases 3 and 4 below. 
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Figures 15, 32, 34, and 36 show the results of CASE 3, in which there 
is no eddy diffusion of the temperature anomaly. It is seen that the 
cold anomaly at the northeast boundary still disappears as in cases 1 and 
2. The western cold anomaly, however, appears to maintain its magnitude 
for a greater duration than in the previous two cases; the greatest dif- 
ference occurring during the 60 to 90 day prediction period. The western 
edge of the warm anomaly moves eastward about the same distance as before 
but its shape is altered and the eastern edge remains unchanged as in 
cases 1 and 2. The major difference observed in CASE 3 is that the maxi- 
mum amplitude center of the warm anomaly only decreases from .6°C ini- 
tially to .4°C at the end of the 90 day prediction. Consequently, the 
results of CASE 3 imply that the anomalous diffusion terms of (39) were 
significant in producing the heavy damping in cases 1 and 2. Since re- 
moving the anomalous diffusion terms of (39) improved the results of the 
prediction considerably, but did not completely eliminate all of the 
damping, CASE 4 was devised in which the diffusion of the anomalous cur- 
rents was also removed. The removal of these additional diffusion terms 
should eliminate some damping because of the dependence of the tempera- 
ture on the currents through the geostrophic relationship as discussed 
previously . 

With respect to dissipation, movement, and amplitude of the anomalies, 
CASE 4 (Figs. 15, 38, 40, and 42) produced results nearly identical to 
those of CASE 3. Therefore, it is concluded that the diffusion terms 
governing anomalous currents in the primitive equations are seemingly 
unimportant in the prediction of sea-surface temperature anomalies. Never- 
theless, CASE 4 represents a more consistent model in that there is no 
diffusion of either the anomalous currents or temperature. 
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If one now compares the results of Namias' prediction (Fig. 4) with 
the results found in this study using a completely dynamical model, it 
is first observed that the 0° isotherm of the warm anomaly has about the 
same configuration and location. Secondly, the location of the center 
of maximum amplitude in the predicted warm anomaly is much closer to the 
observed 90 day position in Namias* study than in this study. However, 
if one closely examines the location of this anomaly center in Figs. 15 
and 16, it is seen that the anomaly center moves eastward almost 15° of 
longitude during the first 30 days and then remains relatively station- 
ary. This observation would imply that the movement of the anomaly's 
center is not due simply to advection alone as indicated in Namias* pre- 
diction. Lastly, the amplitude of the predicted warm anomaly decreases 
in all cases of this study but the amplitude of the same anomaly in 
Namias' prediction remains constant as expected when employing horizontal 
advection alone. Therefore, the vertical advection in this model plays 
an important role in the evolution of such an anomaly. Whether that 
role in the present prediction or any future predictions is simulated 
accurately by the model depends critically on the vertical profile 
assumed for the anomaly. 
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IV. CONCLUSIONS 



This study has shown that the method of prediction using model gener- 
ated climatology along with observed temperature anomalies in a fully 
dynamical numerical model of the ocean circulation can be successful. 

Also, the simple geostrophic initial conditions for the anomalous cur- 
rents is successful away from the equator. However, the quality of the 
above results depends on the precise initial data available. 

In view of this study, the need for temperature data (normal and 
anomalous) in the mixed layers of the ocean is concluded to be the dominat- 
ing factor influencing the prediction of the sea-surface temperature 
anomalies. This appears quite evident from our observation discussed 
above in reference to Figs. 15 and 16. Since this rapid movement of 
the anomaly center cannot be explained by horizontal advection alone, 
it is speculated that a sub-surface warm anomaly may have existed down- 
stream which was brought to the surface. Thus, vertical advection may 
have played a significant role; however, the lack of data beneath the 
surface makes these speculations virtually impossible to investigate. 

The present model, with no salinity, cannot accurately predict the 
evolution of high latitude cold anomalies which do not extend to the 
ocean floor. Therefore, initialization should include salinity; however, 
at present, large scale anomalous salinity data is not available. 

When comparing the results of cases 3 and 4, it is observed that 
the removal of the diffusion terms governing the anomalous currents had 
little effect, if any, on the 90 day prediction. Thus, further 
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predictions using a dynamical model should, for consistency, omit all of 
the anomalous diffusion terms (CASE 4). 

Although the results of this study did not produce predictions which 
corresponded exactly with the observed anomaly field, the use of a fully 
dynamical model in making these predictions has given us a greater in- 
sight as to the role played by vertical advection and diffusion. It is 
therefore believed that further studies should be conducted using this 
fully dynamical model with more accurate initial data and with salinity 
included. 

Additionally, this study uses only climatological data to define the 
atmospheric forcing. However, an anomalous atmosphere would produce 
anomalous heating and wind stresses which could produce anomalous features 
in the ocean. For example, Namias (1965) investigated the effects of an 
anomalous wind stress and showed that this anomalous atmospheric forcing 
may be important in some cases. Thus, further investigations should be 
conducted to determine the response of the ocean to various anomalous 
atmospheric forcing . 
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Igure k. Initial (a) and predicted (b,c,cnd d) sea-surface temperature anom- 
alies after 30, 60, and 90 days. Hatched aroa3 show anomalies ranging between 
-IF and -2F; stippled areas indicate anomalies >1F. Heavy vertical line nark* 
170W. Fron Namias(1972a.fie;.9). 
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Figure 5. Domin of predictions 
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Figure 6. Initial normal currents generated by 11-year integration 
of model. 
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Figure 7. Initial normal temperatures generated by 11-year integration 
of model. 
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Figure 8. 30-day normal currents 
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Figure 9. 30-day normal temperature field 
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Figure 11. 60-day normal temperature field 
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Figure 12. 90-day normal currents 
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Figure 13. 90-day normal temperature field 
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Figure 14. Initial observed anomalous currents. Fall 1971. 
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Figure 15. Initial observed anomalous temperature field. Fall 1971. 
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Figure 16. 30-day observed anomalous temperature field. 
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Figure 17. 60-day observed anomalous temperature field. 
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Figure 18. 90-day observed anomalous temperature field 
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Figure 19. CASE 1: 30-day predicted anomalous currents. 
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Figure 20. CASE 1: 30-day predicted anomalous temperature field 
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Figure 21. CASE 1: 60-day predicted anomalous currents 
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Figure 22. CASE 1: 60-day predicted anomalous temperature field. 
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Figure 23. CASE 1: 90-day predicted anomalous currents. 
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Figure 24. CASE 1: 90-day predicted anomalous temperature field. 
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Figure 25. CASE 2: 30-day predicted anomalous currents. 
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Figure 26. CASE 2: 30-day predicted anomalous temperature field. 
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Figure 27. CASE 2: 60-day predicted anomalous currents 
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Figure 28. CASE 2: 60-day predicted anomalous temperature field. 
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Figure 29. CASE 2: 90-day predicted anomalous currents. 
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Figure 31. CASE 3: 30-day predicted anomalous currents. 
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Figure 32. CASE 3: 30-day predicted anomalous temperature field 
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Figure 33. CASE 3: 60-day predicted anomalous currents 
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Figure 34, CASE 3: 60-day predicted anomalous temperature field. 
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Figure 35. CASE 3: 90-day predicted anomalous currents 
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Figure 36. CASE 3: 90-day predicted anomalous temperature field. 
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Figure 37. CASE 4: 30-day predicted anomalous currents 
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Figure 38. CASE 4: 30-day predicted anomalous temperature field. 
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Figure 39. CASE 4: 60-day predicted anomalous currents. 
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Figure 40. CASE 4: 60-day predicted anomalous temperature field. 
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Figure 41. CASE 4: 90-day predicted anomalous currents. 
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Figure 42. CASE 4: 90-day predicted anomalous temperature field. 
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